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Abstract 

We study properties of the momentum space Triple Pomeron Vertex in perturbative 
QCD. Particular attention is given to the collinear limit where transverse momenta 
on one side of the vertex are much larger than on the other side. We also comment 
on the kernels in nonlinear evolution equations. 

1 Introduction 

The Triple Pomeron Vertex (TPV) in perturbative QCD [1, 2, 3] has attracted significant 
attention in recent years. It is derived from the 2-^4 transition vertex in QCD reggeon 
field theory which represents the high energy description of QCD. In recent years partic- 
ular interest has come from studies of nonlinear evolution equations, e.g. the Balitsky- 
Kovchegov equation [3, 4, 5], where the nonlinearity is given by the TPV. More recently, 
also generalizations of the nonlinear evolution have been considered [6] which contain 
pomeron loops [7]. Again, the TPV plays a central role in these investigations. Whereas 
in many studies and applications it is convenient to use the coordinate representation, it 
is important to understand the structure also in momentum space. 

In this paper we will investigate some aspects of the TPV, starting from the momentum 
space representation of the 2 — 4 gluon transition vertex, from which the TPV vertex has 



originally been derived [1, 2]. Many of the studies of the nonlinear evolution equations 
have been done in the context of deep inelastic scattering where a virtual photon scatters 
off a single nucleon or off a nucleus. In both cases the momentum scale of the photon 
is much larger that the typical scale of the hadron or nucleus, i.e. one is dealing with 
asymmetric momentum configurations. As a first step of investigating the TPV, therefore, 
we will focus on investigating the limit where the transverse momenta are strongly ordered. 
We also review and discuss the nonlinear evolution equation that have been proposed in 
the literature [8, 9], and we comment on the use of a twist-expansion in the low-x limit. 

The paper is organized as follows. In the following section 2 we define the setup 
of our calculation, and we construct the elastic amplitude for photon-photon scattering 
with exchange of a four gluon BKP state. In the large- iVc limit, this reduces to photon- 
photon scattering with the exchange of a pomeron loop. In the remaining part of section 
2 we define the Mellin transform of a pomeron loop and specify our use of 'collinear' 
and 'anticollinear limits'. In section 3 we study the collinear limit of the Triple Pomeron 
Vertex in the large Nc limit. Section 4 contains results of the analysis of the anticollinear 
limit in the large Nc limit. In section 5 we extend the analysis to finite Nc. In section 6 
we derive a hierarchy of nonlinear evolution equations which describe the interaction of 
a photon with a hadronic target. We also show that, in the mean field approximation, 
we obtain a nonlinear evolution equation for the unintegrated gluon density. Section 7 
contains a few comments on the relation of this equation with other nonlinear evolution 
equations described in the literature. We end the paper with a few conclusions. 

2 The 2^4 gluon transition vertex 

The LO momentum space expression for the 2 — > 4 gluon transition vertex has been de- 
rived in connection with the diffractive dissociation of the virtual photon in deep inelastic 
electron proton scattering [2]. More precisely, the process 7* + g {qq + n gluons) + q 




(a) (6) 

Figure 1: Examples of diagrams that contribute to the 2 — 4 gluon transition vertex 
(wavy vertical lines represent reggeized gluons): (a) real emission, (b) a disconnected 
contribution. 

has been investigated in the triple Regge limit. The resulting vertex consists of three 
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pieces (we follow the notation of [11]): 



m - 1 



(J"^'^^(J"^"*^(ACi,K2,ki,k2,k3,k4) 



+<^"^"^(^"^"^V^(/«1, /«2, ki, k3, k2, k4) + /«2, ki, k4, k2, k3) 



where K^i + k;2 = ki + k2 + k3 + k4 = q, and the subscripts a^, Oj refer to the color degrees 
of freedom of the reggeized gluons. It is convenient to express the 'basic vertex function' 
V{ni, K.2, ki, k2, ks, k4) in terms of another function G(k.i, k.2, ki, k2, ks): 



V{ki, K2, ki, k2, k3, k4) = -g* 



K2, ki, k2 + k3, k4) + G(Avi, K2, k2, ki + k3, k4) 
+G{ki, K2, ki, k2 + k4, k3) + K2, k2, ki + k4, k3) - K2, ki + ks, k3, k4) 

-G{Ki, K2, ki + k2, k4, kg) - G{k,i, K2, ki, k2, k3 + k4) - K2, k2, ki, kg + k4) 

+G{k,i, k,2, ki + k2, -, k3 + k4) 



(2) 



This function G(k.i, K2, ki, k2, ks) [13, 16] generahzes the G function introduced in [2] to 
the non-forward direction. This function can again be split into two pieces: 

G{ki, K2, ki, k2, kg) = Gi{ki, K2, ki, k2, kg) + G2{k,i, K2, ki, k2, k3), (3) 

where the first part contains the 'connected contributions' (also: 'real contributions'): 

k2K.f K.2 



Gi{Ki, K2, ki, k2, kg) = + 



(/Ci-ki)2 (K2-k3)2 (Ki - ki)2(,.2 - k3)2 

-(kl + k2 + k3)^ 



(4) 



and the second one takes care of the disconnected ('virtual') pieces: 

2 2 

/G'2(/ci,K2,ki,k2,k3) = -^([a;(k2)-a;(k2 + k3)]5(')('*i-ki) 

+ [cu(k2) - W(ki + k2)]5(2)(Avi - ki - k2)). 

Here a;(k) denotes the trajectory function.: 



(5) 



+ (k-l)2(k-l)^ 



(6) 



The vertex (1) is completely symmetric under the permutation of the four gluons. It is 
infrared finite, it has been shown to be invariant under Mobius transformations [12], and 
it vanishes when or kj goes to zero. 



3 



Qi-q 




Qi-q 



Qo + q 



(a) 




\> + q 



Figure 2: Contributions to the elastic scattering of two virtual photons which contain the 
2 —>■ A gluon vertex (dark blobs represent Green's functions of reggeized gluons): (a) the 
four gluon BKP state; (b) a Pomeron loop. 



This vertex can be used to construct, in reggeon field theory, the selfenergy S, of 
the BFKL Green's function (Fig. 2). In the lowest order contribution to S, we have a 
BKP state between two 2 — > 4 vertices, which contains all pairwise interactions of four 
reggeized t-channel gluons. Its Green's function satisfies the following evolution equation: 



{lj - u;(ki) - u;(k2) - u;(k3) - u;(k4))^< 



(4) {«,:},{<} 



({k.},{k'J) 



Hi) ^'^i 



({k.}{k:}). 



(7) 



where we have used the shorthand notation {kj} = (ki, k2, ks, k4) etc. The sum extends 
over all pairs {ij) of gluons, the kernel i^^2^2^^*^ includes the color tensor faib^cfajuj 



biC 



T^{a,}^{b,} _ 2 f e 
-'^2^2 y J biaicj ca2b2 



(ki - k)2 (ki - k) 



(8) 



and the convolution symbol ® stands for J-^r^- The inhomogeneous term has the form 



(27r)3 

s^'\j:^.-j:K)gi^^'i'^^}i<\{k,},{K}) = (2^)9]^ ^--<^^^Hk.-kD 

1 



(9) 
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Figure 3: Quark quark scattering in the high energy limit of QCD (color singlet exchange): 
(a) two loop correction in the ladder approximation, (b) diagrams with two Triple Pomeron 
Vertices (grey bloobs). 




Figure 4: Different cuts. 

Let us briefly explain how we obtain the correct normalization factors. We consider the 
scattering amplitude of quark quark scattering (Fig. 3b) which, in the center, contains the 
insertion of four reggeized t-channel gluons. We begin with the diffractive cut (Fig.4a): 
the one loop amplitude on the Ihs of the cutting line will be assumed to have an even 
signature exchange and, hence, is equal (up to factor i) to its energy disconinuity. As a 
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result, we start from Fig. 3b with three cutting hues: all horizontal lines are on mass- 
shell. To distribute the color and phase space factors we proceed as follows: using the 
on-shell conditions, we can perform 8 out of the 10 longitudinal Sudakov integrations; 
the two remaining longitudinal variables denote the rapidity of the two produced gluons 
in the central region which, for the moment, we keep fixed. For each closed loop, we 
are left with an integral J (Pk/{27r)^. For each BFKL rung in Fig.3a we have the kernel 
from (8), for each 2^4 vertex in Fig.Sb we have the 2^4 vertex (1) (divided by the 
additional factor l/\/2). Having in mind that our discussion should be applicable also 
to more general diagrams, we retain, for the moment, the general color structure in (8) 
and in (1). When deriving, from Fig.3b, the diffractive cut in Fig. 4a, we write, for each 2 
gluon exchange on the Ihs and on the rhs of the cutting line a statistic factor 1/2!; from 
the compensating factor 4 we absorb V2 into each of the 2 — > 4 vertices. As a result, we 
have, in addition to all other color and phase space factors, the statistic factors 2/(2!)^. 
Invoking now the AGK rules [14, 15], applied to the exchange of four (odd signature) 
reggeized gluons, the other contributions in Figs. 4b and c, the statistic factors become: 

4 

(10) 



2!2! 3! J 4! 

Finally, we use our result for quark-quark scattering and return to the process of our 
interest, photon-photon scattering. Replacing the quark impact factors by photon impact 
factors, inserting BFKL rungs above and below the 2 — > 4 vertices, and inserting pairwise 
interactions between the four gluon lines in the center, we arrive at: 

A{s,t) = 2isTT [ dYs [ dY2 [ dYiS{Y -Y1-Y2-Y3) 
Jo Jo Jo 

7^7^*^ ' ' q - «i, q 

(27r)'^ (^TxY 

•if /n (1^) (2-)^^^^nEk. - q) /n (H) ^^^n\Y.K - q) 

,ya'ia'-a,a,a,a,^^^^ q - A^i; ki, k^, kg, \.^)Q^''^-^^^^'^^''^\Y2, {^iMK'i) 

y (27r)3 (27r)3 K2, K3, K4, q 

.gi2Wlb'^m^y^. q)/!^^..', q - k'), (H) 

where the minus sign in the third line indicates that the four gluon insertion into the two 
gluon Green's function represents a negative correction to the simple ladder amplitude. 
Here s is the squared center of mass energy, Y= ln(s/so) is the total rapidity, Yi, Y2, 
Ys are the rapidity intervals as depicted in Fig. 2, 0"'i"2 denotes the impact factor of 
the virtual photon, ^^^)"i"2'"i"2 jg ^^le BFKL Green's function which satisfies the BFKL 
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integral equation: 



(a; - a;(kO - u;{k,))g^J^ --^^^^{k,}, {k'J) = 

^(2)0a,a.;.i6.(|k,},{k^}) + T^K^^^'^ ® Q^J^ '^^'^^'^^"^iKji^}), (12) 

with an inhomogeneous term analogous to (9). The statistics factor ^ reflects the symme- 
try of the expression under the interchange of the four gluons. In eq.(ll), the selfenergy is 
defined by lines 3-5, i.e. the convolution of the two 2 — > 4 vertices with the BKP Green's 
function between them. As a convenient simplification, we approximate the four gluon 
state by two noninteracting color singlet ladders (Fig. 2b): this configuration represents 
a 'pomeron loop'. It is easy to find the combinatorial factor of a system where two pairs 
of gluons form bound states. We have three possibilities of pairing two gluons to form 
bound states out of four gluons. This yields the factor 1/2!. In this configuration we have 
a pomeron loop topology. The result reads: 

A{s, t) = 2^S7^ [ dY^ I dY^ [ dYi S{Y -Y1-Y2- Y^) 
Jo Jo Jo 

■ q - /.)^(^)»i«-«"«^'(y3; K„ q) 

4 / (0^ / ^^^"'^'"^'^"^"^"^"^(-' ^ - - r' + r + q) 

■ / ^^^^^)^'^""''''^^^' ^^' ^'^' ^)iPSf^^'^'''''iy2; ka, k^, r + q) 

.gi2W(W,^Y^■, /.', q)/i^H/.', q - /.'), (13) 



where 



Xaia2Xbib2 

501026162 _ (^d\ 



is the color singlet projector. These projectors act on the color tensors of the 2 — > 4 
vertices, turning the pairs of color labels (0102), (&1&2), (03O4), (bsb^) into color singlets. 
Comparison with (1) shows that this projection operator, when acting on the first term, 
leads to a factor 1, whereas the remaining terms come with the weight factor j^2_i ' in 
comparison with the first term, they are color suppressed. This large- TVc approximation 
turns the 2 — > 4 vertices into the Triple Pomeron Vertices (TVP). 

In the following we shall focus on the pomeron loop (13) and investigate, for zero 
total momentum transfer, q = 0, the kinematic limit where the momentum scale of the 
upper photon is much larger than the lower one. This implies that, at the upper TPV, 
the momentum from above, Ki, is larger than the momenta from below, ki, ka, and the 
loop momentum r ('coUinear limit'). Conversely, for the lower TPV we have the opposite 
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/m(7) 

' > 


— 1 




^ -Re (7) 



Figure 5: Singularities in the 7 plane. 



situation: the momenta k'^, kg, and r are larger than k'^^ ('anticolhnear hmit'). Let us 
become a bit more formal. We expand the amplitude of Fig. 2 in powers of Qq/QI 
('twist expansion'). The object of our interest is the self-energy of the Pomeron Green's 
function, T,{ki,k[). In Eq.(13), T,{ki,k[) is defined to represent the lines 3-5, i.e. the 
convolution of the two TPV's with the two BFKL Green functions between them. It has 
the dimension k^, and it is convenient to define the dimensionless object 



with the Mellin transform: 

poo 

s(7) = / det{e){ey-\ (15) 

Jo 

The inverse Mellin transform reads: 

„2^ _ f '^(l.2^-^^ 



S(n = y^^(^)""^W' (16) 

where k = and the contour crosses the real axis between —1 and (see Fig. 5). Our 

analysis will then reduce to the study of the singularities of the function 2(7). The twist 
expansion corresponds to the analysis of the poles located to the left of the contour in the 
7 plane: the pole at 7 = — 1 is the leading twist pole, the pole at 7 = — 2 belongs to twist 
4, and so on. As we have already said before, for the upper TPV in Fig. 2, the analysis 
of this twist expansion requires the 'coUinear limit', for the lower TVP the 'anticolhnear' 
one. 
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3 The collinear limit 



In this section we are going to study the coUinear hmit of the TPV. The ordering of the 
transverse momenta is the following: |k;i| = |k|^|ki|, |k2|, |k3|, |k4|. We therefore expand 

in powers of |ki|/|k|, |k2|/|k|, |k3|/|k|, |k4|/|k|. In our investigations we will be interested 
in attaching color singlet objects to the vertex, and we project (1) onto the color singlets. 
In the limit Nr ^ oo we obtain: 



paia2bib2 pa3a4b3b4'^a[a'2;aia2a3a4 k" k^ k2 k3 k4) = 

V{1, 2, 3, 4) + ^ {V{1, 3, 2, 4) + ^(1, 4, 2, 3)) 
where V{1, 2, 3, 4) = ^(k, — k; ki, k2, k3, k4). The first term will be denoted by 



vlo'^^''^(l, 2, 3, 4) = (5«i.«^5^'i.fe5^'3,b4 



TT 



the second and third ones by 



3, 2, 4) = s<A^b,,b2^b3M 



\/27r 



- 1)^ 



V^(l,2,3,4), 



rni,3,2,4) 



(17) 



(18) 



(19) 



etc. 



3.1 The real part 

Let us begin the analysis by expanding the real part of the G function (3) in the collinear 
limit. As we are going to limit ourselves to the forward case we use the simplified notations 
G(k, -k; ki, -ki - kajka) = ^(kijks) and \/(k, -k; ki, k2, ks, k4) = l^(ki, k2, ks, ki). 
With these definitions the Gi function (4) reads: 

k^k^ k^k^ rk, + ko^^i^"^ 



(k-ki)2(k + k3 



In the collinear limit the momenta of the outgoing gluons satisfy the conditions |ki|«|k|, 
|k3|«|k|. Performing the expansion in |kj|/|k| up to fourth order terms we obtain: 



G'i(ki,k3) = 2k2 



kik3 ki-k 



kvk 



(k^ + 2kik3) + -V (k? + 2kik3) + 



+ 



k^ .2krk 2 



)M (k2 + 2kik3) + 



,2k,k. 



k^ 



k^ ^ k^ 

+(ki + kg) 



(k^ + 2kik3) + (21) 



22k ki 2k k3 
k^ 



+ ... 



The first term is the twist-two contribution: 



G'i(ki,k3)^=' = -2k 



kik3 



(22) 
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With analogous expressions for the other Gi functions in eq.(2) we find that the sum 
of all twist-two pieces vanishes. The next two terms on the rhs of eq.(21) vanish after 
averaging over the azimuthal angle of k. Finally we are left with the twist-four piece. 
After averaging over the direction of k we find: 



G'i(ki,k3)^=" = 2k 



k' 



(23) 



Prom eq.(2) we obtain for the twist-4 piece of the real part of the TPV: 

vi:;2!'\k,MMMr='= 

xa\ .g; ;c6i.&2 xbzM 9^ 01.2 • k2 kg • k4 + ki • U ka • kg + ki • kg k^ ■ k^) 

^ ^ ^ k^ ' ^ ' 

where the superscript r stands for the real emission. This expression is the master formula 
for the twist-four contribution. 

In the next step we are going to attach BFKL ladders to the pairs of gluons (ki, k2) 
and (k3,k4). Since the presence of a momentum transfer across the BFKL ladder would 
cause loss of a logarithmic contribution, we limit ourselves to the forward directions: 

ki = -k2, k3 = -k4. (25) 

Putting ki = 1, k2 = —1, ka = m, k4 = — m we obtain: 

VIS">(1. -1. m. -m)- = ^.i.»i^>..>^^»..>.-|^^2k- ^P"-";;r (26) 

Now we multiply the vertex by propagators for the lower gluon lines and convolute with 
2 — 2 transition kernels (cq. (8)). Our goal is to find, from the convolutions of the vertex 
with propagators and kernels, the maximal number of logarithms. To do that wc should 
act on the twist four contribution of the vertex with the twist four evolution operator, 
which, in our case, is the product of two BFKL kernels in the twist-two approximation. 
Let us compute the coUinear approximation to the BFKL kernels which, when convoluted 
with the TPV, will give a logarithmic integral. The expression for the emission part of 
the BFKL kernel is: 



^(qi,q2;ki,k2) = -N^g 



(ki + k,f "i'^' 



(k2 - ci^y (ki - Qi)^ 



(27) 



The factor — A^c replaces the color tensors in eq.(8), since we have projected on the color 
singlet state. Assuming zero momentum transfer and qf ^kg we get: 

K^g''N,2kl (28) 

where, in order to simplify the notation, we have skipped the arguments of the kernel K. 
Using this approximation in the formula for Fig. 6a, we get the following expression for 
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(a) 



Figure 6: (a) The TP V with two BFKL interactions attached to it. (b) The TP V with 
gluon ladders. 

the convolution of the vertex with two BFKL kernels (one kernel for each two-gluon pair 
below the vertex): 

f'^' dH r""' d^m 2kl 2kl 4(2(1- m)2 - l^m^) 
M (2^Ag (2^1^^ k^ 
= 0, (29) 

where kg is the lowest momentum scale which we do not specify at present. We notice 
that the integrals of I and m are logarithmic, but - what is most striking - the angular 
integral over the angle between I and m renders the Triple Pomeron Vertex to vanish. 
It is straightforward to iterate the convolution with BFKL kernels (Fig. 6b), and, as a 
result, we arrive at the conclusion that - after averaging over the azimuthal angles - the 
twist-four part of the TPV vertex gives zero contribution. 



3.2 The virtual part 

So far we have investigated contributions coming from the real part of the vertex. What 
remains are the disconnected parts. In order to investigate logarithmic contributions of 
the virtual pieces we have to convolute them with an impact factor at the upper end of 
Fig. lb. To deal with infrared finite quantities it is convenient to work with the impact 
factor of the photon. The function G2(ki,k2 + k4,k3) = G'2(ki,k3) (5) in the forward 
direction reads: 

k3) ^ -k'gL (,„ _ k) + >.>^J^^'^'(k3 + k,) . (30) 

The photon impact factor (for transversely polarized photons) has the form [20, 21] 

0.,.,(k,Q) = 6<<a.a.^Y.< ^^^P ^'' ' '^T^ " "J^^ ^ (31) 

„ Jo p(l -p)g2 + r(l -r)k 
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where and denote the (negative) photon and gluon virtuahties, resp. In this expres- 
sion p denotes the longitudinal component of the quark loop momentum (in the Sudakov 
decomposition), while the second integration variable r is a Feynman parameter. For our 
investigations we are interested in the twist expansion. To perform the twist expansion of 
the impact factor one has to perform the Mellin transform with respect to k^/Q^. With 
the Mellin transform (15) one finds: 



0aia^(k, g) = y ^^\^) M-'M^ 



and we obtain: 



r(3 + 7) r(i + 7) r(-7) r(-7 + 2) 
1 + 7 -7 r(-7 + |) 



1 + 7 



(32) 



(33) 



Turning to Fig. 2b, we are interested in the following ordering of momenta: ^ |k| » 
^ iQol- To analyse the twist-4 term of this kinematic region we need, for the upper 
impact factor, the twist-4 term. Closing, in eq.(32), the contour of the 7 integration to 
the left we obtain the following collinear expansion of the photon impact factor: 



14 4, 

¥-3^H^ 



where C = ejasCKem- As it is well-known, the twist- four term has no logarithmic 
enhancement. 

For later purposes wc also list the results for the lower impact factor: we close the 
contour to the right and find: 



0<.lai(K',go) = i"'"^(«',Qo) = S'''"^C- 
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1(01 

5 U'' 



(35) 



For our twist-4 analysis of Fig. 2b we will need the second term on the rhs. 

Returning to the upper impact factor and concentrating on the twist-4 piece, we now 
easily see, by simply counting powers of momenta, that the virtual contributions of the 
TPV cannot contribute to the maximal power of logarithms. Namely, beginning with 
the impact factor above the TVP, we have the power k^ which cancels the two gluon 
propagators attached to the impact factor. From the G2 functions we find another power, 
k^, which, through the 5-functions, turns into m^, 1"^, or (m ± I)''. When dividing the 
region of integration into the two parts m -C 1 and 1 -C m, the terms with (m ±1)^ turn 
into 1^ or m^. Below the TPV we have the pairs of propagators, 1/m'^, and 1/1'^, and 
there is no m (or 1)- dependent contribution from the BFKL kernels. Combining these 
momentum factors, we therefore obtain only integrals of the form 



1^ 
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i.e.none of the integrals is logarithmic (this argument remains unaffected if we include the 
logarithms from the G2 functions). Hence, within the leading-log approximation, also the 
virtual part of the TPV is zero. 

Let us emphasize that our search for the 'maximal power of transverse logs' is exactly 
what is required for a consistent twist-four analysis. In order to obtain this maximal 
power (i.e. one power for each transverse momentum loop integral), we had to start, 
at the upper impact factor with the twist-four term. Including a BFKL-ladder between 
the impact factor and the TPV forces us to take also the twist-four approximation of the 
kernel, i.e. instead of the leading twist approximation in (27), terms of the order 0{\<?fy). 

2|2 

Next, at the TPV we searched for terms of the order ^^if-, and, finally, for the two BFKL 
kernels below the TPV, again the twist-2 approximation (27). It is only this sequence of 
approximations which provides one logarithm for each loop, i.e. otherwise we loose one 
(or more) powers of logarithmic enhancements. Our result then says that one coefficient 
in this sequence of terms, namely the TPV, vanishes and thus makes the twist-four term 
in the twist expansion (in the leading-log approximation) disappear. 



3.3 Generalization to all higher twists 

The main result of the previous subsections - the absence of coUinear logarithms in the 

case of angular averaged BFKL ladders - can be generalized to all orders of powers of 
l/Q^- We return to the function V in (2) which is expressed in terms of the functions Gi 
and G21 and we average over the angles of m and 1. First Gi. 

n I ^_ k'l' k^m^ k^(l + m)^ 

^) - (k^ip + (k + m)2 " (k-I)2(k + m)2- ^^^^ 

Let us denote the first term in this formula by A, the second one by B and the third one 
by C, the angle between I and m by a, and the angle between m and k by /3 (the angle 
between 1 and k then equals 27r — a — /3). For the integrals over a and /3 we find: 



1 Z"^'" Z'^'' k^l^ 



" " (2 

and: 

To compute the integral over C we spht C = Ci + C2 + C3 into three pieces. The 
corresponding integrals are: 



27r ]^4j2 

dadj3-—2 : 2\ 2 (^9) 



(27r)2 Jo '■^ (k^ - 21kcos(a + (3) + f){\<? + 2mkcos /3 + m2) ' 
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= J_ riR 2k'|l||m|cosa 



1 f'^'^ k^in^ 
^'''^j2^Jo - 21kcos(a + (3) + \'){k' + 2mkcos/3 + m2) ' ^^^^ 

The results of the integration are: 

l^k^ 

= 112 ,2|| 2 1.2r (42) 

|1 — k ||m^ — k I 
-Sl^m^k^ 

- |j2 _ ^2||j^2 _ k2|^i2 ^ ^ _ ]^2|^(^2 + 1,2 ^ |^2 _ k2|)' (^S) 



m^k^ 



Il2 1 2| I 2 1 2r (44) 

1 1 — k 1 1 m"^ — k I 

The total contribution is given by summing up I a, Ib, Ici, Ic2 ) ^c^ The result can greatly 
be simplified if we consider special situations. For instance, if k^ -C 1^, m^, we may drop 
the absolute value signs. Adding all terms we obtain: 

V- , 2l2mV - 2lV - 2bV + 2k6 ,2 

In all other cases: k^ l^,m^, m^::§>k^3>l^, or 1^3> k^:»m^ the sum of all terms gives 
zero. Therefore the final result can be simply written as: 

/ dad(3G,{\,rn)^2\^e{\^-\^)e{xa^-\^), (46) 

where the factor l/(27r)^ comes from averaging. 

Let us now perform the angular averaging of the disconnected pieces of the ^(l, m) 
function. We have: 

G.(W, ^ (in (I^^'^'d - k) + ln^^<^'(m - k)) . (47) 

To compute the integral over angles we split the region of integration. In the case when 
|m|^|l| the first term gives 



ID 



2iT j2 j2 

dad(5 In ^ = In ^ , (48) 



(27r)2 Jo 1^ + + 21m cos a 



whereas the second one vanishes. In the case when |I|^|m| we obtain zero from the first 
term and 



1 m2 

= - — — / daddln -5 = In ^5- (49) 

{2ny Jo l' + m2 + 2ImcosQ; ^ ^ 
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from the second one. We combine the two cases in the following way: 



^2(1, m) 



87r2 



ln-^^(m2- 12)^(2) (1 



m 



TTI 

k) + ln^^(l^ 



m )5(')(m-k) ) . (50) 



Putting all pieces together and including the remaining G functions we arrive at the 
angular-averaged form of V: 



1 



(27r) 



/■27r 

/ dad(5V{\i, -k; 1, -1, m, -m) = 4 

Jo 



2k'e(l' - \^)e{ra^ - k") 



(51) 



The presence of the ^^-functions forbids all collinear configurations, i.e. there is no 
expansion in inverse powers of k^. The physical meaning of this result is the following: if 
the two pomerons entering the vertex from below have smaller momenta than the pomeron 
from above, they cannot resolve it and cannot merge because they do not feel 'color' and 
the vertex vanishes^. In the language of a twist-expansion, our result states that, in 
the leading-log approximation, not only twist-four, but all higher twist terms are zero, 
provided we restrict ourselves to the large-A^c hmit, and we use only the BFKL ladders 
with conformal spin zero below the TPV. 

4 The anticollinear limit 



4.1 Real part 

Let us now investigate the anticoUinear limit of the 2 — > 4 vertex (Fig. 7). In contrast to 
the collinear limit where the BFKL ladders below the vertex had to be in the forward 
direction, the anticollinear configuration allows for a nonzero momentum transfer across 
the BFKL ladders above the vertex. The momentum transfer here, as we will see, does 
not lead to a loss of a logarithm. We arc interested in the limit |w|«|wi|, |w2|, |w3|, |w4|. 
To study the real emission part of the TPV it is convenient to rewrite the Gi function in 
the form: 



Gi(Wi, W2 -|- W3, W4) = 

2 (Wi + W4 



+ 



(1_2^ + ^) (1 + 2^ + ^) 

|2 1 



-W 



wfw4 



fl _ + + 2^^ + 



(52) 



similar result has first been noticed in [25]: however, the disconnected pieces have been missed. 
The result (51) agrees with the form given in [2]. 
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Here wc have used the momentum conservation Xlj'^j — 0- The expansion parameters 
are |w|/|wi| and |w|/|w3|. Performing the expansion we obtain up to second order: 



-2^ - 2^ - 2w 

W4 



G'i(wi, W2 + W3, W4) = 

2 Wi ■ W4 



2 + 2 



^ W-Wi 

wf 



W- Wi 



^ W-W4 

' 2 



W-W4 



(53) 



Using (4) we obtain for the leading term of the TPV: 

,r{a'} 



V 



LON, 



(p, -p - r, q, -q + rY^^^ding ^ ^a'^a'^^WM^hzU J ^ g^2w'. 

i\„ — 1 2 



(54) 



One easily sees that this term does not provide logarithms in the momentum scale. The 
subsequent terms in (53) vanish after averaging over the angle of w. Therefore, in order 
to get, after convolution with BFKL kernels in the subsystems (12) and (34), the required 
logarithmic contribution we need to consider, in (53), terms of higher order. After av- 
eraging over the angle of w, and after summing, in (4), over all the Gi functions, the 
resulting contribution is the following: 



^LOJVe(Wl,W2,W3,W4) 



V2: 



TT 



-g w 



(J«'l'«2(j''l'''2^&3,fc4 

"'^-'■^ "'IW3 W2W3 

Wi-(Wi+W2) W2-(Wi + W2) W3-(W3+W4) 



Wi-W3 W2-W3 Wi-W4 W2-W4 



W?Wo 



WnWo 



WIW4 



WoWa 



W4 ■ (W3 + W4) 

wf(Wi + W2)^ w|(Wi+W2)2 W|(W3 + W4)2 W4(W3 + W4)2 

(Wi + W2)-(W3 + W4' 



(55) 



(Wi + W2)2(W3 + W4)2 

To proceed further we need the anticoUinear limit of the BFKL kernel. Using (27), 
setting ki = Wi, k2 = W2 and requiring that |qi|^|wi|, |w2| we obtain: 



K = - g^Nc2wi-W2. 



(56) 



As already mentioned before we are interested in the maximal power of logarithms in the 
momentum scale; this leads to the particular momentum configuration, where the mo- 
mentum transfer across the BFKL Pomerons in the subsystems (12) and (34) is nonzero, 
(note that below the vertex we are still in the forward direction). We set: wi = p, 
W2= — p — r,W3 = q, W4 = — q-|-r. In order to obtain, after convoluting with BFKL ker- 
nels in the subsystems (12) and (34), logarithmic contributions, we have to consider the 
following momentum-ordered configurations: 

• The configuration where |r|^ |p| ^ |q|. BFKL kernels and propagators are of the 
form: 

A^c/2p-r, -A^c/2qT (57) 
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and 




resp. In order to render all transverse momemtum integrations (in p, q, and in r) 
logarithmic, we need, from the TVP, terms proportional to they are obtained 
from the first term in (55): 

^LON, IP, -p-r,q, -q + rj - -0 1 ^ at2 _t9 ^ -rZo l^yj 

C Mr ^ 

Combining these expressions and performing the integrals we obtain: 

(K. A-,) » VIS:'"-- = ^«'.4i'.'.i'.'._|lLA,|_|!_^to ^) ' . (60) 

where |wo| is the momentum at the upper end of the BFKL kernels, specified by the 
condition that it should be smaller than the momentum scales |1| and |m| which 
were considered in the coUinear limit of the upper TPV. Convoluting this expression 
with the impact factor </){fe/}(w) = 5'''i''2|C^ InQo/w^ below the vertex yields: 

iK^ K,) ® Vl^^^f^'—' ® = 5'^''^5'-^''^V2.Nf-C^^ (in ^) ' . (61) 

Here, in order to get the logarithmic contribution, we took, in (35), the next-to- 
leading term in the anticollinear expansion of . In the configuration |r| ^ |q|^|p| 
the same result is obtained. 

• Repeating a similar analysis in the case when |q|^ |p|^|r|, we find, from the last 
term on the rhs of (55): 

The same contribution is obtained from the region |p|^ |q|^|r|. 
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• Finally, there are the regions |q| ^|r| ^|p| and |p| ^|r| ^ |q|. For the first case we 
use the first term in the second line of (55) and obtain: 

(K, K,) ® ® ha'} = 6'-'^^6'-^^''^V2nNf-C^^ (in ^) ' . (63) 

The second region gives the same contribution. 

4.2 Virtual parts 

Let us now analyze the contribution coming from the virtual parts of the vertex in the 
anticoUinear limit. Again, we are looking for the maximal power of logarithms. We begin 
with the region |r| ^ |p|, |q|. Using (2), (30) for the virtual parts of the TPV, (57) for the 
BFKL kernel, (58) for the propagators above the vertex, and (35) for the lower impact 
factor we immediately see that none of the G2 functions allows for three logarithmic 
integrals. The same observation holds for the regions |q| ^|r| ^|p| and |p| ^|r| 3> |q|. 
We are then left with the region |p|, |q| ^ |r|. From the BFKL kernels and from the 
propagators wc find the denominators 1/p^ ■ 1/q^, and we therefore need the factor 1/r^ 
from the propagators below the TPV. They can come only from the first term in the 
function ^2(1 + 2, 3, 4): 

4/2 2 \ 

/G2(l + 2, 3, 4) = a,^ ( In <5(w + W3 + W4) + In ^ <5(w + W4) ) , (64) 

27r V (W3 + W4)2 w| J 

and from analogous terms in 6^2(1 + 2, 3, 4), 6^2(1, 2, 3 + 4), 6*2(2, 1,3 + 4). Convoluting 
these G2 functions with the BFKL kernels and with the impact factor, and setting the 
lowest momentum scale equal to Ql, we find: 

(i^i K2) VS'^^> = -5'-'^^5'-'^V2nNf-C^^ (in . (65) 

Let us shortly summarize our results for the large-A^^c limit, before wc continue the 
finite Nc analysis. Our goal was to find those terms of the twist expansion of the TPV 
which, after convolution with the BFKL kernel, would generate the maximal possible 
power of transverse momentum logarithms. In the coUinear case (upper TPV), we had 
to restrict the BFKL ladders below the TPV to the forward direction, and we therefore 
expected to find, from the m and 1 integrations, two logarithms. After the convolution 
with the upper impact factor, a third logarithm should appear. What we found is that 
the coefficient of this maximal number of logarithms vanishes, both for the connected 
and for the disconnected parts of the TPV. In the anticoUinear case we had to include 
the integral over the momentum transfer across the first BFKL kernel. After convoluting 
these integrals with the lower impact factor, we expect four logarithms. In fact, we found 
these logarithmic contributions, both in the real and in the virtual part of the TPV, and 
they came from different regions of ordered transverse momenta. 
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This completes our twist-4 analysis of the one-loop Pomeron self-energy of the BFKL 
Pomeron (Fig. 2b). We have found that the upper TPV vanishes at the twist-4 point, 
whereas the lower one provides nonzero contributions. 



5 Finite 

5.1 The coUinear limit 

In this section we are going to investigate contributions to the vertex in (17) that are 
suppressed in the large Nc limit. Repeating our analysis of the previous sections we 
obtain for the first subleading piece: 

yr{a'},{b}f^ o o 4^^=^ _ 

_ xa'.a'oxbub2xb3M .4 ,,2 4 ("kj • kgka • k4 + ki • kakg • k4 + ki • kik2 ■ ks) 

^ ^ (iv2-i)2^ ^ ^4 ■ m 

Substituting ki = 1, k2 = —1, k^ = m, k4 = — m we obtain: 

V 

(7V2_i)2^" k 



V:Sf ^(1, m, -1, -m)-- = ri"^5^-^M^3,..V27r ^4^211!^. (67) 



The convolution with the two BFKL kernels gives: 

VslT'^-' ® (^1 K.) - ^"'"'^"'''^"''-^^^.^^^ (l- 1) ' • (68) 

Convolution with the impact factor gives: 

« v^^f^-' « (K. = |l (,„ I)' . m 

The same result holds for the second subleading part. 

For the virtual corrections to the TVP the situation is the same as for the leading- A^^ 
part: by simply inspecting the powers of transverse momenta, we find that the integra- 
tions over m and 1 are not logarithmic, i.e. they cannot generate the maximal power of 
logarithms. 

5.2 The anticolUnear Umit 

Here our starting expression for the real part of the TPV can be taken directly from the 
rhs of eq.(55), by interchanging W2 and W3. For the first nonleading piece we have: 

ViF^(Wi,W3,W2,W4)— ^ = 
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Wi-W2 W2-W3 Wi-W4 W3-W4 

2 2~ 2 2~ 2 2~ 2 2~ 

W^W2 W2W3 W|W| W3WI 



Wi-(Wi+W3) W3-(Wi+W3) W2-(W2 + W4) W4 ■ (W2 + W4) 



(70) 



W^(Wi+W3)2 W§(Wi + W3)2 W|(W2 + W4)2 W4(W2 + W4)2 

(Wi + W3)-(W2 + W4) 
(Wi + W3)2(W2 + W4)2 

The analysis is analogous to the leading- A^c case. In detail we find: 

• For |r|^|p|^|q| the logarithmic contribution comes, on the rhs of eq.(70), from the 
second term of the second line. We obtain: 

(^1 ^2) ® Vi2f ^ ® <Pia'} = ^'-''^'-''^.Nclc^^ (in ^) ' . (71) 

The same result holds for the region |r| » |q| ^ |p|, taking in eq.(70) the first term 
of the second line. This result is same as in (61), except for the suppression by 
- 1. 

• |q|^|r| ^ |p|: here we use the first term on the rhs of eq.(70) and obtain 

{K^ K,) V:i2'}'^ ct>ia', = ^''^'^^''''^.Nllc^^ (in ^) ' . (72) 

The region |p| ^ |r| 3> |q| gives the same result. It coincides with (63), but is 
suppressed by A*"^ — 1. 

The regions |q|:»|p|:»|r| and |p|:»|q|:»|r| do not contribute to the maximal number of 
logarithms. 

Finally we come to the virtual parts of the A^^c-suppressed parts of the TPV. Repeating 
the analysis, carried out for the virtual part of the leading- A'^c piece, we find no contribution 
to the maximal number of logarithms. The final result for the anticoUinaer limit of the 
A'c-suppressed part of the TPV, therefore, is given by the real piece alone. 



6 Nonlinear evolution equations 
6.1 General evolution equations 

Let us now make some use of the TPV in QCD rcggcon field theory. To be definite let 
us consider deep inelastic scattering on a hadronic target (a single proton or a nucleus). 
We define color singlet i-channel states of n reggeized gluons {n even) in the Heisenberg 
picture which are labeled by color and momentum degrees of freedom: 

\n) = -^al{\^,)...aiSKm 
\/n\ 

= |ki, ..k„;ai, ..an). (73) 
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The normalization is: 

[a„(k), ai,(k')] = (27r)=^k^5(k - k')5„„, (74) 

and 

= J] n ((2^)'<^(1^ - k^k.^^-^"^^) , (75) 

cr{n) i=l 

where the sum extends over the permutations of outgoing gluons. The unity operator is 
given by: 

oo n /" (/•^k- 1 

^Y)){ri = ^Y[ — ^^^i,..k„;ai,..an)(ki,..k„;ai,..a„|, (76) 

n n 1=1 ' i 

where the summation on the left hand-side includes also the integration over the contin- 
uous degrees of freedom. 

We assume that the target state, at some initial rapidity, can be written as a super- 
position: 

oo 

\P)-J2^nil). (77) 

n=l 

The rapidity evolution of this (color singlet) state is given by: 

e'"\p) = \p{y)), (78) 
The Hamiltonian consists of several pieces(Figs. 8-10): 

H — i?2-+2 + -f^2-»4 + -f^4-»2 + -f^2-»6 + -f^6-*2 + ■■■ ■ (79) 

The first term denotes the case where, inside the n gluon state, only one pair of gluons 
interacts, in the second term one pair splits into four gluons etc. The matrix elements of 
H2-,2 are expressed in terms of the BFKL Hamiltonian: 

n' . . 

{n'\H2-.2\n) = 5r,n' J2 U<cfcaf^a, \ i^2^2(k„ k, ; k',, ^){27rfS{K + k, - k^ - k;.) 

i>j=l L 

+ (a;(k,) + u{k,)) k2k|5(k, - k[)6{k, - k;.)| + (k: ^ k;, a: ^ a'^) 

n' 

J2 n(2^)'^(k^-k;)^a.a; (80) 



(n'-2)\ ^ 

where the second term in the square bracket stands for the symmetrization of the outgoing 
two-gluon state. This kernel corresponds to the BKP interaction in the color singlet 
state. All the other terms are presently known only for the special case where not only 
the total n gluon system but also the interacting subsystem belongs to the color singlet 
representation. 
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{n'\H2^2\n) = 




Figure 8: Matrix element given by eq.(80). 



{n'\H,^2\n) = 



{n'\H2^4\n) 




Figure 10: Matrix element given by eq.(82). 



In particular, the second term contains the 2 — > 4 transition vertex: 



n n 



V^y^^-'^-'^'iK, k^, k;, k;; k„ k^ + (k, ^ k^, a, ^ at) 



{n'\H2^i\n) = 5nn'+2 ^ ^ 

s>t=l i>j>t>i=l 

1 " 

(27rf<5(k, + k, + k, + k,-k,-k,)^--— J2 n (27rmkp-k;)<5'^''<, 

<T(n-4) pj^i,j,l,r 

(81) 

whereas the third term allows for four gluon to fuse into two gluons: 



{n'\H4^2\n) = Snn'-2 ^ 5^ 



ya.ar,a,a,a,a^^^>^^ ^. ^ ^ k^, a', ^ 



(27rf5(k, + k,+k, + k,-k;-k;)— -— 5] n (27rmkp-k;)5'^^"^- 

a(n'-4)pi^i,j,l,r 

(82) 

The next two terms on the rhs of (79) belong to the Pomeron two Odderon vertex 
[22] (restricted to the color singlet channel) and to its inverse, resp. They will not be 
discussed further. Higher order kernels (indicated by the dots) have not been computed 
yet. Let us define the n reggeon wave function component of the target at rapidity y in 
the following way: 

V'^>(y,ki,k2...k„) = (n|e^^|p). (83) 
Upon differentiation with respect to y we obtain: 

^ n' 

= Y,{n\Hyi')i^lf'\ (84) 

n' 

This defines an infinite set of coupled equations. It cannot be closed because, for instance, 
the equation for the two gluon wave function receives contributions coming from the four 
gluon wave function: 

= {2\H2-.2\2)^Pr' - {2\H,^2mr' (85) 
(the term proportional to (81) vanishes since it requires zero gluons in the initial state). 



6.2 The nonlinear equation for the unintegrated gluon density 

In order to reach further simplification, we take the large- A^c limit. In practice this implies 
that we group the n gluons into n/2 color singlet pairs (Pomerons) and associate with each 
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pair a color singlet projector: this projector acts on the color tensors of the interaction 
Hamiltonians and leads to color weight factors of the interaction kernels. In particular, 
in the 2 — > 2 Hamiltonian the color tensor faia[cfca'.aj is replaced by the color factor —Nc, 
and in the 2 ^ 4 Hamiltonian, i72^4, the 2 ^ 4 vertex V"""*'"*"^""'"" reduces to the 
function V^(k'j, k^, k'^, k^) (cf. eq.(l)). The evolution equations have to be reformulated 
in terms of N states of gluon pairs: each pair carries two momentum variables, q and k: 
q denotes the total tranverse momentum of the two gluon state, and k, q — k are the 
momenta of the two constituent gluons. The state consisting oi N — n/2 such pairs is 
defined as: 

\N) = -^At(q,,kO...At(q^,k^)|0) 

= |(qi,ki),..(q^,k;v)) (86) 

(we use capital letters to distinguish the pair-basis from the reggeon basis) . The normal- 
ization follows from 

Kq,k),At(q',k')] = (27r)V(q-k)2%-q')5(k-k'), (87) 

in analogy with the reggeon states. We write the Hamiltonian as 

H = Hi^i + H,^2 + H2-.i, (88) 

where 

{l\H,_,i\l) = 7Vc|i^2-.2(ki, k,-; k^, k;)(27r)35(k, + k,- - k^ - k^) 

+ (^(k,) + ^(k,)) K^k|5(k, - k:)5(k, - k;)| + (k: ^ k;.) (89) 



and 



(l|^i-2|2) 



V{K kj, k^; k„ kt) + {K ^ k*, ^ at) 

^3 



(90) 



{27ry6(ki + kj + k; + k^ - k^ - k^). (91) 



The amplitudes in this basis of gluon pairs are defined in analogy with (83). 

Next we invoke the mean field approximation and make the following factorizing 
ansatz: 

^2{y, ki, qi - ki, k2, qa - ka) = ^i{y, ki, q^ - k,)^,{y, ka, qg - ka). (92) 

This ansatz can be justified for a large nuclear target. It allows to obtain a closed equation 
for *i: 

= - -^(l|i/2-.i|2)*i*i. (93) 
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To obtain the BK equation for the unintegrated gluon density let us define the off-diagonal 
unintegrated gluon density via: 



T{yMM) = MvMM) = (l|e-^^|p). (94) 
Using (80), (81), (89), (93), (94) we obtain the nonlinear evolution equation: 



dlnl/x J (27r)3 ' ' l'(q-l)2 

,2 dH d'^m , , , q q N 

^ '(2^ (2^^^*^' -k + q; 1, -1 - 2 + r, m, -m - - - r) 

T{,v. I + r, 1) ^(,r. f - r. m) 

X " 



(95) 



Y{-\ + I + r)2 m2(-m + | - r)2 ' 

The momenta entering the vertex from below are labeled by k^=l, kg^— 1— q/2-|-r, kg^m, 
k4=— m — q/2 — r. The variable r stands for the loop momentum. In [11] it has been 
shown that this equation coincides with the Balitsky-Kovchegov equation, provided the 
solutions JF belong to the Mobius class of functions (i.e. the Fourier transform vanishes 
when the two coordinates become identical). We make the assumption that the coupling 
to the proton goes via the form factor (with momentum transfer r) 



e 



nr,i?) = ^^, (96) 

(where R has the meaning of the proton radius), and for J-{x, r, k) we make the ansatz: 

J'ix, r, k) = T{x, k)F(r, R). (97) 

Then the integration over r on the rhs of (95) will be restricted to small values < 
Now we restrict ourselves to zero momentum transfer, q = 0, which corresponds to the 
integration over the impact parameter, and, as further approximation, we put r = at 
the TPV: in the dipole language, this means that the typical dipole size is assumed to be 
much smaller than the impact parameter b. This allows to carry out the r integral, and 
one easily sees that the function J-{x, k) satisfies the somewhat simpler equation: 

dJ^{x,k) f dH , , ^ ,:F{x,r) 



1 f d'l d'm J^{, 
-"2^7 (2^(2^^^'^'-'^'^'-^'"^'-"^)^ 



J (2.)3^^1' 1'^' ^) 1^ ^^^^ 



d\nl/x 

jF(x, 1) J^{x, m) 



In the next step we perform the integrations over the azimuthal angles of 1, m, and k. 
Denoting the integrated function J^{x, k) by f{x, k^): 

fix, k') = ^ y k) (99) 
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with 



xg{x, k^) 



/ 



J2 



(100) 



and using our result (51) for the angular averaged TPV, the nonlinear equation reads[8, 9]: 

2 




1 

2 



Ik^-l^l 



/(^,n 



(101) 



When applying this equation to the scattering of a virtual photon on a nucleus we 
return to the question raised at the end of the introduction, the question of the most 
dominant gluon configurations. In the DGLAP approach one has strong ordering in 
momentum, i.e virtualities of gluons closer to photon are larger than those closer to the 
target. In the nonlinear evolution equation one then would expect that, at the kernel of 
the nonlinear term, the upper momenta, k, should be larger than the lower ones, 1^ and 
m^. However, making use of our results for the coUinear limit of the TVP and of the 
structure of the angular averaged vertex, we arrive at the somewhat surprising conclusion 
that the momenta are ordered in the opposite direction. In more physical terms, the 
recombination of two smaller gluons ends up in a larger gluon. This suppression of softer 
gluons below the nonlinear term may explain why, in numerical solutions of the angular 
averaged BK equation for the unintegratcd gluon [9], the BFKL diffusion into the infrared 
region is absent. 



7 Comparison with other equations 

As we have mentioned before, the nonlinear equation (95) coincides with the Balitsky- 
Kovchegov equation. In [11], the Fourier transform of (95) has been computed, and it 
has been shown that, in the class of Mobius functions, it agrees with the BK equation. 

Alternatively, one can start [8, 9, 10] from the Balitsky-Kovchegov equation for 
the dipole scattering amplitude in coordinate space, and compute the Fourier transform 
to momentum space. The connection between the momentum space gluon distribution 
J-'{x, q, k) and the dipole scattering amplitude is: 



>(a;,q,k)-^^^^^k(k ^) J J e e 



(102) 



where Xqi^Xq — Xi, and b= (xi + X2)/2 is the impact parameter. Our steps of approxi- 
mation described after (95) are equivalent to the factorization ansatz 

Ar(xoi, b, x) = N{^oi,x)S{h) (103) 
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and to the assumption that, in the Bahtsky-Kovchegov equation, all dipole sizes are 
much smaller than the impact parameter b. With these approximations one arrives, after 
angular averaging and integration over impact parameter b, at the nonlinear equation 
(101). 

Returning, once more, to the issue of the twist expansion, we have to conclude that the 
nonlinear BK equation, when restricting to solutions with conformal spin n — 0, receives 
all its contributions from 'anticoUinear' terms. This, in connection with corrections to the 
single- ladder approximation at small x, makes the usefulness of a twist expansion some 
what doubtful. 

Let us finally comment on other versions of nonlinear evolution equations. The first 
nonlinear evolution equation which was a milestone in physics of saturation is the Gribov- 
Levin-Ryskin, Mueller-Qiu (GLR-MQ) equation [25], [26] (eqn(2.41) in [25], and eqn.(30) 
in [26]), obtained in the double- logarithmic approximation: 



d'^xg{x, k^) cusNc 
ain(l/x)aink2 ~ TT 



xg{xy)-C:^^[xg{xX)r 



(104) 



(the constant C is not the same in the two papers; however, for our discussion this is 
not essential) . This equation can be rewritten in terms of the unintegrated gluon density 

fixyy. 



df{x,k') N,as 
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fix, I') - C 



a: 



k2 



(105) 



dhil/x 

The linear term coincides with the BFKL kernel in the collinear approximation. The 
nonlinear term should be interpreted as the TPV at the collinear limit. Its physical 
interpretation would support the strong ordering (collinear) picture discussed at the end 
of the previous section: momenta above (/c^) are larger than below [k'"^) the nonlinear 
interaction. Our analysis, however, does not agree with this form of the nonlinear term. 
The structure of integrals is totally different. In particular, we have come to the conclusion 
that, after angular averaging, the TPV does not contribute to the collinear limit. 

The GLR paper [25] also presents another nonlinear equation (eq.(2.108)), derived 
from summing up, at small x, single logs of the fan diagrams. It is written directly for 
the unintegrated gluon density which, in the GLR notation, differs from our definition 
(eq.(lOO): 



xg{x, 



fc2 



cik''$(x, k'^) 



(106) 



This equation is an attempt to generahze the BFKL equation to the physics of dense 
systems, and its form is quite close to our equation (101): 



d^xX) 
dh\l/x 



°° d\^ 
T 



ll^-k^l 



1 



-Qtpv 



( 



a. 



AtvR^ V47r 



(107) 
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where grpv is the local approximation of the following TPV vertex: 



(108) 

This vertex contains the same 6'-functions as in (51), and it thus supports the physical 
picture describes at the end of the previous section. On the other hand, the detailed 
analytic form of the vertex is different from (101); in particular, it does not contain the 
disconnected pieces which, in the original derivation of the 2^4 vertex, can be traced 
back to the reggeization of the gluon (there are also differences in the prefactors) . 

Despite these differences in the detailed form of the nonlinear equations it may very 
well be that, as far as the gross features of saturation are concerned, the qualitative 
behavior of solutions will be similar. It would be interesting to study this in more detail. 



8 Conclusions 

In this paper we have investigated the momentum space triple Pomcron vertex. In partic- 
ular, we have studied its collinear and anticollinear limits. This question arises naturally if 
one studies nonlinear corrections to the linear BFKL evolution in deep inelastic scattering 
at small x: one expects that, at least on the average, transverse momenta decrease when 
moving from the photon to the proton. In a first step one is then led to consider the limit 
of strong ordering. Restricting ouselves to solutions with conformal spin zero, we have 
shown, for the simplest example of a fan diagram with one triple Pomeron vertex in the 
large- iVc limit, that there is no contribution from the configuration of strongly ordered 
gluons. Beyond the large- iVc limit such contributions exist. 

We have also constructed a set of evolution equations for the interaction of a photon 
with a nuclear target, which, in the mean field approximation, reduces to a nonlinear 
evolution equation for the skewed unintegrated gluon density which, in the forward region, 
agrees with equation obtained in [8, 9]. We have also compared our momentum space 
expression for the nonlinear evolution kernel with different other versions discussed in the 
literature. We agree with the BK equation, but we find disagreement with other earlier 
versions of nonlinear evolution equations. 

Interpret at ing our results in terms of twist, we have shown that the BK-equation, 
when restricted to solutions with conformal spin zero, receives all its contributions from 
'anticollinear' configurations, quite in contrast to the expected ordering of transverse 
momenta. 

We also hope that our analysis will help to analyse further the contributions of pomeron 
loops. 
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